SISSA 38/97/EP 
NSF-ITP-97-034 
February 1, 2008 



A^-species Stochastic models with boundaries 

and 

quadratic algebras 



F. C. ALCARAZQ 

Departamento de Fisica 
Universidade Federral de Sao Carlos 
13565-905 Sdo Carlos, SP, Brazil 

S. DASMAHAPATRAg 

Department of Mathematics 
City University 
London ECl V OHB, UK 

V. RITTENBERG0 

SISSA, Via Beirut 2-4 
Trieste 1-34014 Italy 
and 

Physikalisches Institut 
Universitdt Bonn 
D- 531 15 Bonn, Germany 



^Partially supported by CNPQ and FAPESP-Brazil 
^Partially supported by the EPSRC grant GRJ 25758 

^Work done under partial support of the EC TMR Programme grant FMXR-CT96- 
0012 



ABSTRACT 



Stationary probability distributions for stochastic processes on linear 
chains with closed or open ends are obtained using the matrix product 
Ansatz. The matrices are representations of some quadratic algebras. The 
algebras and the types of representations considered depend on the boundary 
conditions. In the language of quantum chains wc obtain the ground state of 
A^-state quantum chains with free boundary conditions or with non-diagonal 
boundary terms at one or both ends. In contrast to problems involving the 
Bethe Ansatz, wc do not have a general framework for arbitrary N which 
when specialized, gives the known results for A'^ = 2; in fact, the N = 2 and 
N > 2 cases appear to be very different. 



1 Introduction 



The aim of this paper is to present in a systematic way the apphcation of 
quadratic algebras to obtain the steady state probabihty distribution of one- 
dimensional stochastic processes with boundaries. In this section we first 
present the problem and describe what is already known about the subject 
giving the relevant references and then proceed by giving the contents of 
the next sections. The list of possible physical applications of our results 
include interface growth Q , boundary induced phase transitions IH, |3|, ^ ^ , 
the dynamics of shocks or traffic flow Q . 

Consider a linear chain with L sites. On each site k we put a discrete 
stochastic variable (3^ taking values from {0, 1, A'^ — 1}. For each link k 
between sites k and + 1 we give transition (intensity) rates T/j*^^^^^ giving 
the probability per unit time for the transition 

Ik Ik+i Pk Pk+i 

o — o - o — o 

These are the bulk rates. At the left end of the chain (site 1) we also 
consider transition rates 

71 01 

o - o 

given by the N x N matrix L^^ . Similarly, at the right end of the chain (site 
L) we take processes 

o - o 

with rates i?^^ . The matrices L and R give the boundary rates. We consider 
Markov processes in continuous time, which implies that the rate matrices 
are intensity matrices, which have the property Q 

N-l N-l N-1 

p/3fe/3fe+i _ _ -pA/3fe+i. r/3 _ _ — _ R/3 

^/3fc/3fe+i~ J- 7fc7fc+i . -^/3 - -^7' ~ ^7 ' 

7fe>7fc+i=0 7=0 7=0 

(7fe,7fc+i)#(/3fc,/3fe+l) l¥=P 7¥=0 

(1.1) 
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Sometimes it will be convenient to write the intensity matrices (for ex- 
ample R) in an alternative form: 



Rjs — -R/37 



(1.2) 



indicating by the column (7) the initial state and by the row (/?) the final 
state. Equations (1.1) just state that the sum of the matrix elements on each 
column of an intensity matrix is zero. The time evolution of the probability 
distribution P{(]i, P2, ■ ■ ■ , Pl', t) ='■ Pi{P}'yt) is given by the master equation 



d_ 

dt 



L-l N-1 



k=l 7fe.7fe+l=0 



(1.3) 

Solving equation (1.3) is equivalent to finding the wave-function of an imag- 
inary time Schrodinger equation (see Q and references therein for the no- 
tations) which is obtained as follows. We consider an orthonormal system 
of states ^ 

(1.4) 



and a basis in the space oi N x N matrices E°'^: 



(1.5) 



By i?^^ we denote the matrix E^^ acting on the /c site. The probability 
distribution P({/3};t) is mapped into a ket state: 

/ P(0,0,..,0;t) \ 

P(o,o,..,i;t) 



IP)) = Y.n{m)\i 

V ^'(A^-l,iV-l,-,A^-l;i) ; 
and the master equation (1.3) implies the Schrodinger equation 

d 



(1.6) 



where 



dt 

L-l 



PI 



-H \P) 



-f^ — E ^k,k+l + Hi + Hl, 



(1.7) 
(1.8) 



fe=i 
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Hk,k+i = -T'^sE'l°'Ek+v (1-9) 

Hi = -L'^E^'' and Hl = -R'^eI''. (1.10) 

Above, and subsequently, the summation is invoked on any pair of repeated 
indices (raised, lowered, or mixed). 

In this paper we are interested in finding the stationary probabihty dis- 
tribution Ps{{P}) of the master equation (1.3), i.e. the ground state wave 
function \Ps)) of the Hamiltonian (1.8): 

H \Ps)) = 0. (1.11) 

Since H is an intensity matrix with positive rates the ground-state energy is 
zero. Unless stated otherwise we consider unnormalized probability distri- 
butions. The interest in knowing the stationary distributions of stochastic 
processes with boundaries is illustrated by the extensive list of references on 
the subject which can be found for example in [^]. If the bulk rates give a 
Hamiltonian with a certain algebra or superalgebra as symmetry (see ref. 



[10| for many examples) the boundary terms might break this symmetry. 
Recently Krebs and Sandow [11 1 have proven the following remark- 



able theorem. With T, L, R matrices as above, take matrices Da [a = 
0, 1, — 1) and A^ matrices Xq, (a = 0, 1, A^ — 1) acting in an auxiliary 
vector space which satisfy the quadratic algebra relations 

r'^^DaD^ = D^Xs-X^Ds, i-f, 6 = 0,1,. ..,N-l). (1.12) 

Consider now a ket state |0) and a bra state (0| in the auxiliary space (we 
use the same notation suggesting the "vacuum" of quantum mechanics for 
reasons which will be apparent later). If these states are chosen such that 
the following conditions are fulfilled: 

{0\{Xa-L^aD/3) = and {X^ + R^^D^M = 0, (1.13) 

then 



PsiiP}) = (oi n^Aio) 

k=l 
L I N-\ 



:i.i4) 



(o| n E I |o) 

k=\ \/ift=0 



^Observe that the quadratic algebra (1.12) is unusual since the X^'s appear lin- 
early. (See refs. and for a discussion of quadratic algebras in the mathematical 
literature.) 
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is a stationary solution of the master equation (1.3). Alternatively, let us 
denote by u^j. (n = 0,1, . . . , N — 1; k = 1, . . . , L) basis vectors in the vector 
space associated with the k^^ site {{u^^),^ = 5^^u for any k) on which the 
matrices act (see eqs. (1.9-10)). The ground state wave-function of the 
Hamiltonian (1.8) has the expression 

L N-l 

\Ps)) = {o\x{{Y,D,^u,^m. (1.15) 

k=l /ifc=0 

Actually, the Krebs-Sandow theorem is even stronger; the matrices T"^, 
and in (1.8) do not have to be intensity matrices. If the Hamiltonian 
has an eigenvalue zero, then the wave function is given by (1.15). More- 
over, Krebs and Sandow have shown that the algebra (1.12) exists, giving a 
representation for the matrices and X^. This representation fulfils the 
second part of eq. (1.13) but not the first. 

In order to compute the ground state wave-function one needs the "vac- 
uum" expectation values of monomials of the form 

mi\Dll...Dl^^\0), (1.16) 

which are obtained from eqns (1.12) and (1.13). In order to find the wave- 
function of the two-site problem one considers "vacuum" expectation values 
of monomials of degree two which are obtained from a system of linear 
equations. The solution of this system is not necessarily unique since the 
ground state might be degenerate. On completing this exercise the L = 3 
"vacuum" expectation values (1.16) are considered, and so on. This entire 
process, however, is more complicated than diagonalizing the Hamiltonian 
by brute force, and there should be a better way of solving the problem. 
This is indeed possible if the intensity matrices appearing in the Hamilto- 
nian satisfy certain conditions. At this point, we shall mention only what 
is known in the N = 2 case. In reference [0] Hinrichsen et al guessed a 
four-dimensional representation for the matrices Dq , Di , Xq and Xi for a 
specific choice of intensity matrices and derived certain concentration pro- 
files. In this specific case the problem can also be solved directly since the 
Hamiltonian can be written in terms of free fermions and it was this solution 
which inspired the guess. This reference is remarkable since it is the first 
place where the relations (1.12-16) appear. 

We shall now describe the cases where the representation theory is un- 
derstood. Consider the case where the only non-zero bulk rates are Fq^ and 
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r^o but leave the surface (boundary) rates L^, Lq and Ri, Rq arbitrary. This 
is the asymmetric exclusion model with open boundaries. One can choose 
Xi = —Xq = 1 (the c-number 1) in equations (1.12) and (1.13) and af- 
ter a linear transformation D ^ Y which includes constants related to the 
boundary rates, bring the relations to the form 

[Yo,Y^]= Yl c^%Ys+ E c% + coi. (1.17) 

7,<5e{0,l} 76{0,1} 

where -B] = AB — BA and c^^ = . The parameters c^^ ^c^ and cqi 
depend on the bulk and surface rates and we also obtain 

(o|yi = o, yo|o) = o, (o|o)/o. (i.is) 



The quadratic algebra with two generators (1-17) is well understood and 
the Fock representations defined by eq. (1.18) are known [Q, ^. It was in the 
work of Derrida et al that Fock representations of quadratic algebras were 
used for the first time to find ground state wave-functions and to compute 
correlation functions. As one can notice from (1.17), (1.18), the "vacuum" 
expectation values of a monomial of degree n in the Y^s is determined by the 
"vacuum" expectation values of a monomial of degree n — 1 and one of degree 
n — 2 in the Y^s. This implies the existence of recurrence relations between 
stationary distribution functions ({/?}) of lattices of different length L 
(which is also the degree of the monomial in question). The first calculation 
of correlation functions in this model were performed using these recurrence 
relations without using the algebraic approach ^ |3| . 

One may try to take more bulk rates than the purely diffusive or hop- 
ping ones just described. The problem is that (1.12) with X^s chosen to 
be c- numbers gives three equations for the two generators Dq and Di, while 
one is enough to determine all "vacuum" expectation values as we have 
just seen. This implies relations between the rates. One trivial case cor- 
responds to one-dimensional representations for the Da (see Q) when the 
(connected) correlation functions vanish. Even two-dimensional representa- 
tions don't exist, since the relations on the rates obtained are incompatible 
with their positivity |jl5|. Another possibility is to not restrict the X^s to 
be c-numbers, but (say) to take 

Xa = Aa^D^ + Xa, (1.19) 

where A^'y is an arbitrary 2x2 matrix and x^s are arbitrary c-numbers. 



Using the arbitrariness of the matrix elements A^-y, Krebs et al |15] have 
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shown that it is possible to include more non-zero rates than just the hopping 
ones (r^o and Fq^) and get a two-dimensional representation of the quadratic 
algebra, and hence obtain non-trivial correlation functions. 

We would like to mention that the algebraic approach can be used not 
only in sequential processes described by the master equation, but also for 
parallel dynamics |jl^, 0, |l^, where the algebraic problems to be solved 
are identical. Also, it is well worth recalling that matrix-product states first 
appeared in the literature for periodic boundary conditions in the work of 



Hakim and Nadal |1£], and have also appeared in various publications by 
Fannes, Nachtergale and Werner ||2^ who also considered Bethe lattices. 
Parallel work was done by the Koln group [^]. The algebraic problem for 
periodic boundary conditions is very different since in this case, the wave- 
function is given in terms of traces of monomials and not in terms of inner 
products in Fock states. Quadratic algebras can be used in this case too 



(see Refs.m, ll,!!). 

In the present paper, we address the problem of stochastic processes 
involving N > 2 species. We shall consider, for the first time, not only 
the case of open boundary conditions, but also closed or mixed boundary 
conditions since the algebraic approach applies to these cases as well. In 
section 2, we first review some properties of intensity matrices which we 
shall use in the following sections. In section 3, we consider the problem of 
the linear chain with closed ends (free boundary conditions). The boundary 
matrices R and L are identically zero and we shall consider solutions where 
Xa = 0. This case is not only interesting on its own, but it represents the 
natural first step before considering the chain with open ends where R and 
L are non-zero. (The choice = is compatible with eq. (1.13) and will 
take us to an interesting class of quadratic algebras.) This implies 

r'^^D^D^ = 0. (1.20) 

We shall call a quadratic algebra which contains only quadratic terms like in 
eq.(1.20) a polynomial algebra. (Obviously the algebra defined by eq.(1.17) 
is not one.) Equation (1.20) describes relations among monomials of degree 
two. The number of relations depends on the rank of the bulk intensity 
matrix and using eq. (1.14), we get the probability distribution for the two- 
sites case. In order to have a solution for more sites, one needs consistency 
conditions which imply that the two ways of relating the cubic monomials 
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DaDpDj and D^DpDa as depicted in the figure below 

/ \ 
D^DpD^i D^DpDc, (1.21) 

\ / 

give the same result. The requirement of commutativity of the above di- 
agram imposes constraints on the bulk rates and solves the L = 3 case. 
In all the examples presented in the paper, we have checked that once the 
relations (1.21) are satisfied, no supplementary conditions on the rates arise 
from quartic (and in some cases) higher degree monomials. Similar condi- 



tions appear in quantum group structures |25, 26, ^ and it is no accident 
that various quantum planes or superplanes are solutions of the consistency 
conditions in some of the cases that we consider. Let us observe that for 
the closed chain case, one does not need to take average values (0| • • • |0) in 
eq. (1.14); the ground state wave- functions (the ground state is often degen- 
erate) and therefore the probablity distributions are the coefficients of the 
independent monomials. In section 3 we also answer the following question: 
if the ground state can be obtained by the present algebraic methods, is the 
Hamiltonian integrable? We shall show that the answer is in general, no. 

In section 4, we consider the case where the left end of the linear chain is 
open (L^ 7^ 0) and the right end is closed {R^ = 0). We again take = 
but now we need to add the condition: 

Li{{)\Dp = Q (a = 0,l,..,iV-l) (1.22) 

to equation (1.20), which comes from eq. (1.13). This implies new consis- 
tency conditions between the bulk and the left-boundary intensity matrices, 
which come from the equations 

L^{0\DpD^ = 0. (1.23) 

The consistency conditions depend on the rank of the L matrix. Note that 
in this case, one has to drop the |0) symbol and retain the (0| symbol in 
equation (1.14), where the stationary probability distribution is expressed. 
Again there are no Fock representations in this case. 

Fock representations appear in section 5, where we consider the linear 
chain with both ends open. This problem is much more difficult than the 
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preceding ones. We therefore limit ourselves to the case = 3 only, since 
this case is complex enough and essentially different from the N = 2 case. 
We also restrict our attention to the case of diffusion only (in the bulk), i.e. 
the only non-zero bulk rates are T^^. We look for representations for the 
case where the Xq,'s are c- numbers. To our knowledge, only one example 
of this kind was known up to now It is possible to perform a linear 
transformation to bring (1.12) and (1.13) into the form 

7,5e{0,l,...,A'-l} 76{0,l,...,Af-l} 

where 

with = 3, and the Fock conditions derived from eq.(1.13). We review the 
properties of the algebra (1.24) and stress that their properties are essentially 
different if > 2. We also show that the Fock conditions for N > 2 are in 
general too numerous. This leads to a careful separation of cases depending 
on the rank of the boundary intensity matrices. 

In section 6, we show that if all the minors of the boundary matrices are 
non-zero, one can have representations of the quadratic algebra of dimen- 
sion at most two with the bulk and boundary rates lying on some algebraic 
variety. In section 7, we consider the case in which only one minor of and 
one of R'^ are non-zero. We give the representations of the algebra in this 
case. Other cases (three or more minors non-zero, for example) are not pre- 
sented since the paper is long enough even without them. In Appendix A, we 
present simple physical processes where the formalism may be applied. (The 
results in the different sections of the paper may be used for several other 
applications.) We consider the problem of spontaneous symmetry break- 
ing in a two species exclusion model with asymmetric diffusion proposed by 
Evans et al |^. More boundary and bulk rates are considered than in the 
original model. This extension might allow a better understanding of the 
physics of the problem. 

In section 8, we consider the case in which all the principal cofactors are 
zero (the boundary intensity matrices have rank one). We are going to give 
four examples for this case. 

The reader who might be tired of this long introduction is invited to go 
directly to section 9, where a guide to finding the new results is presented 
along with a number of open questions. 
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2 Some properties of intensity matrices 



An N X N matrix M whose elements My are such that the sum of all its 
entries on each column vanishes, i.e. 

N-l 

J2Mij = 0, {j = 0,1,..., N-l), (2.1) 

i=0 

is called an I-matrix. I-matrices are closed under multiplication and form 
an algebra. 

For stochastic problems, one considers I-matriccs which have to satisfy 
the further restriction that the off-diagonal elements are real and non neg- 
ative, since they are interpreted as probability rates of certain processes. 
Such an N X N matrix M with real entries is called an intensity matrix 
if 

1. all its off-diagonal entries are non- negative, Mij >0,ij^j, 

2. and the diagonal elements are negative with 

Mu = -J2Mij. (2.2) 

For any matrix A with entries Aij, if the Laplace row expansion of a 
determinant is written as dctASik = J2j ^ij^kj where Aij is the cofactor 
corresponding to the element of A. From this we can infer an impor- 

tant property of any zero-column sum matrices; namely their cofactors for 
each column are equal in magnitude, and in the case of intensity matrices, 
all the cofactors have the same sign. Let us denote the column-j cofactor 
of an intensity matrix M by the corresponding symbol in calligraphic font 

My 

It therefore follows that a system of linear equations 

N-l 

MmnXn = (2.3) 

n=0 

has a solution given by 

Xn = ^Mn (2.4) 

for any constant (independent of n) ^ if the rank of M is iV — 1. 
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Consider two intensity matrices F and G, and their sum 



E = F + G, (2.5) 

which is also an intensity matrix. The following identity is true ioi N = 1,2 
and 3 only: 

N-1 N-1 

^ ^ (-^nm Gjirnj^m ~ (-^ 1) ^ ^ {-^nmGm Gnrri'^ m)-i (^•6) 

m=0 m=0 

where as noted earlier, the calligraphic iS, Q with subscript m denotes the 
cofactors of F, G for the m*'' column. The last identity is useful for the 3 
species case. Since we are going to use them often in the next few sections, 
we give the expressions for the three distinct cofactors of a 3 x 3 intensity 
matrix F with matrix elements F^: 

:Fo = F^Ff + F^F^+F^Fl 

Ti = F^F^ + F^F? + F^F^, (2.7) 
= F^F^ + F^F^ + FlFl 



3 Steady states for a linear chain with closed ends 
(Ground states for quantum spin chains with 
free boundary conditions) 

As mentioned in the introduction, we consider ground state wave-functions 
of the form 

L N-1 

= n ( E (3-1) 

k=l tJ-k=0 

where the matrices satisfy equations (1.20) and the consistency condi- 
tions (1.21). The type of wave-functions we get depends on the rank of the 
bulk intensity matrix T"^. Let us first consider the case where the rank is 
the maximum possible, i.e. iV^ — 1, and all principal minors are non-zero. 
We denote by Qai3 the cofactors of T^g (see section 2). Using equation (2.4) 
we get 

DaDp = ^D^Ds. (3.2) 
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There are two ways in which the cubic monomial DaDfjD^ can be related 
to D^DpD^: 



(i) D^DpD, = ^D,D,D, = ^^D,D,D, 

y^iv ypa 



(3.3) 



Comparing (i) and (ii) in equation (3.3) (and setting p = v) we get: 

QapQu-y = GavGp-y (3.4) 

and hence the condition on the cofactors (on setting ly = a and /3 = 7) is 

Ga/3 = ^\jQaaQpi3 = G/3a- (3.5) 

The +(— ) is taken if N is odd (even). This implies that the D^s can be 
taken as c-numbers: 

Da = CVG^, (3.6) 

where ^ is a constant. The ground state is unique and the (connected) 
correlation functions manifestly zero. 

The whole structure of the possible solutions for the ground state wave- 
function depends on the rank of the bulk intensity matrix We shall 
restrict ourselves to the cases N = 2,3. Before proceeding further, we would 
like to consider some questions about integrability. We are going to find 
that for some intensity matrices we can find ground-state wave-functions. 
In general the ground is degenerate, which might imply the existence of 
some symmetry in the Hamiltonian. What is this symmetry? On the other 
hand, one can also ask the (possibly related) question, if the ground states 
are obtained from polynomial algebras, is the Hamiltonian integrable? In 
general, the answer is negative. Exact integrability requires more than the 
exact expression for the ground-state wavefunction. 

In certain cases, the defined by 

e/c := (3.7) 

where Hk^k+i is the Hamiltonian density given in equation (1.9), satisfy the 
relations of a Hecke algebra: 

efcCfeiiefc — Ck = ek±iekek±i - ek±i 

[ek.ej] = 0, for|fe-j|>2 (3.8) 
4 = (9 + g~'^)efe. 
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In this case, one can find a spectral parameter dependent solution R{u) of 
the Yang-Baxter equation via " Baxterization" |2^. One can check that 
the generators satisfy further relations among themselves which define 



various quotients of the Hecke algebra [30|. For our purposes, it is sufficient 
to state that the i/^^fc+i belong to such quotients that are classified by a 
pair of natural numbers (P, M) . For each quotient labelled by (P, M) the 
Hamiltonian whose density is given by the appropriate defines a spin 
chain |l[ with UqSu{P\M) (which is the Schur-Weyl dual of the (P,M) 
Hecke quotient) as its symmetry algebra. (For further details, see jlOt] .) 
A more general criterion for integrability at the level of Hamiltonian 



densities was introduced by Reshetikhin |32] for the case of i?-matrices whose 
dependence on the spectral parameter is of difference type. By expanding 
the R and row transfer matrices in powers of the spectral parameter, it was 
shown to be necessary that the Hamiltonian densities satisfy the following 
relation 

[ek + efc+i, [efc, e^+i]] = Wk - Wk+i, (3.9) 

where Wk is a tensor product of identity matrices at all sites on the chain 
and an arbitrary matrix at the /c*^ site. 

• N = 2 

We now consider the N = 2 case. We shall on occasion refer to the 
state "1" as that denoting the presence of a particle and the state "0" to 
its absence. If all the minors are non-zero, and the intensity matrix is of 
maximal possible rank (i.e. 3) equation (3.6) applies. We consider the cases 
where the rank is less than 3. Let us first consider the case 

roo _ pOO _ tOO _ tOI _ pio _ t11 _ n /-q 1n^ 

^ 10 — -"^ 01 — ^ 11 — ^ 00 — 00 — 00 — {^.LU) 

where the bulk intensity matrix effectively reduces (from a4x4)toa3x3 
matrix with rows and columns labelled by (1, 0), (0, 1) and (1, 1). Let Q'l g, 
G'q i and G'l i be the cofactors of this new intensity matrix. In order to have 
a polynomial algebra we must have the condition 

0[,o = 0'o,i, (3.11) 

and the algebra is 

D^Do = DoDi = ^Df. (3.12) 

The ground state for L sites is now doubly degenerate. One state corre- 
sponds to Dq and the other contains all the monomials which have at least 
one Di, and can be brought by eq.(3.12) to be proportional to Df". 
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Next we consider the situation where we have only diffusion of particles 
and their annihilation if two particles are on neighbouring sites. This is 
attained by choosing the following rates to vanish (the rest being non-zero) : 

pOO _ -pOO _ tOO _ -pOi _ -pio _ -pOi _ pio _ n a 

^ 10 — 01 — 11 — ^ 00 — 00 — 11 — ^ 11 — I^J.IOJ 

In this case the algebra is 

DoDi = q^DiDo; Df = 0, (3.14) 

where 



/If 



(3.15) 





-nOO 
10 





plO 
00 





\ 












1-11 

^ 01 






T^OO 

^ 10 





rio 
J- 00 







V 





Til 





J- 01 


/ 



If we choose the time scale such that rQ^r^g = 1, the set of Hk,k+i sat- 
isfies the relations for the (1, 1) quotient of the Hecke algebras centralizing 
Uq{su{l\l)). The ground state for L sites is doubly degenerate corresponding 
to the words ^^i*^ DiDq~^. 

Yet another example of an intensity matrix which gives rise to a quadratic 
algebra is 



(3.16) 



which gives the following quadratic algebra: 

1 pOO pOl 

DoDi = ^Dj and D,Do = q^Dl where = = (3.17) 

9 ^ 00 ^ 01 

after solving for the associativity constraints coming from the cubic terms. 

Yet another possibility is to set more rates to 0, in addition to those 
which led to equation (3.14). Upon setting FgJ, Tq\ and T\q to zero as well, 
which means we are only considering diffusion, we are led to 

pio 

DoDi = :^DiDo. (3.18) 

^ 10 

The ground state is now (L + 1) times degenerate - each wave-function 
has a given number of Dis and the coefficients of the tensor products of 
the u^j^-'^s in eq.(3.1) correspond to the g-deformation of the symmetrizer 
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corresponding to the Young diagram with one row and L boxes 
^ ~ y^rQ^/r^Q. The quadratic relation 

DqDi = (3.19) 

corresponds not only to the case where F^q is the only non-zero rate, but 
more generally is non-zero for all a, f3. Again the ground state is (L + 1) 
times degenerate. Since our intention is to generalize the problem to > 2, 
it is useful to present visually the non-zero matrix elements which correspond 
to the different quadratic relations. The rows and columns below are labelled 
(0,0), (0, 1), (1,0) and (1,1) consecutively. 



/ 


□ 


□ 


□ 






□ 


□ 


□ 


□ 




□ 


□ 


□ 


□ 


V 


□ 


□ 


□ 





eq.(3.6) 



V • 



□ 


□ 


□ 


□ 


□ 


□ 


□ 


□ 





eq.(3.12) 



□ 
□ 



□ □ 

□ □ 
• □ 



eq.(3.14) 



□ 
□ 



□ 
□ 



V 



eq.(3.18) 





■ □ 




\ 




□ • 


□ 




□ 


• □ 






V • 


□ • 


□ 


J 


/• 


□ • 








□ • 








□ • 






V • 


□ • 


• ) 





eq.(3.17) 



eq.(3.19) 



We could not obtain algebras with non-scalar Dq and Di for any other 
choice of locations in the intensity matrix for the non-zero rates than those 
depicted above. 

• N = 3 

We now consider the = 3 case. We are not going to list all the possible 
quadratic relations for the two-site problem which are compatible with the 
3-site problem (i.e. the cubic relations are compatible with the associative 
application of the quadratic relations) but list a few example which do. We 
shall, in some cases, interpret the state variables 1, 2, as those of two species 
of particles and holes respectively. First, we consider the case in which the 
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two-site bulk intensity matrix decouples into three intensity matrices (the 
other rates being zero): 



-nOl 
10 
■nOl 
10 



i 01 
~^ 01 



r02 

20 

r02 

20 



and 



-nOO 
J- 12 



1-12 
00 

^12 



roo 

21 



r20 

02 

r2o 

J- 02 



p21 
00 
t21 

12 



(3.20) 



pl2 

21 



12^ 



(3.21) 



We denote by J^qq, J\2 and J2\ the cofactors of the 3x3 intensity matrix. 
The quadratic relations obtained are 

D,,Di=^DiDq D2Do = 6DoD2, 
DiD2 = [3Dl aD2Di=[5Dl. 

The cubic relations obtained by considering words in the -D's of degree 3 are 
compatible with the quadratic relations provided 

J21 

There are (2L + 1) wave-functions with energy eigenvalue zero. They are 
given in terms of the monomials 



a 



7 



/? = — ; 
^00 



7 



rio r02 

^ 01 _ -'^ 20 (o 

roi ~ r20- i^o.zz; 
^10 -'^02 



.0 'D^, D^-^D^. (3.23) 

This gives an explanation as to why it was possible to calculate the proba- 
bility distribution exactly in [^3| for the particular case where interchange 
of particles is forbidden (Fgf = = 0)- Interpreting the matrix entries 
in eq. (3.20) as transition rates involving two species of particles 1 and 2 
we notice that the difference in the number of particles of types 1 and 2 is 
conserved under this stochastic dynamics. 

We now consider the intensity matrix of the block diagonal form 





■nOl 
00 


r02 
00 
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-nOO 
01 


-nOl 
J- 01 


1-02 
01 






















■nOO 
02 


■nOl 
J- 02 


r02 
02 































plO 

-■^ 10 


pii 
10 


t12 
^ 10 






















plO 

-■^ 11 


pii 
11 


t12 
^ 11 






















plO 
12 


pii 
J- 12 


yl2 
12 































1-20 
20 


1-21 
20 


1-22 
20 






















r20 
21 


1-21 
21 


1-22 
21 




I 

















7-20 
J- 22 


p21 
J- 22 


1-22 
22 


/ 



(3.24) 
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and we denote the cofactors of the 3x3 block containing TH by G'^-^ , where 

a = 0, 1,2 labels the {a + 1)*'^ column within each block. The algebra is 
defined by the six relations 



D^Df, = ^Dl, (3.25) 



where 



f a 
J a 

/l = £k. = ^ /O ^ ^(2) ^ ^(0) 



(3.26) 



f G"^ G"^ 

T~G^~ C^' ^(0)^(1)^(2) - ^(o)t'(i)t'(2)- l-^-^O 

■/l 5^(1) '^(2) 

The ground state is 3-fold degenerate, the independent words being Dq, 
a,ndD^. This de generacy arises since the type of particle located at the first 
site (0, 1,2) is unchanged by the stochastic dynamics. 
An interesting example with 4 quadratic relations is 

T%DoDi = Tl\DiDo 

r20n n — r02n n i^o.zoj 
Dl = 0. 

These relations are obtained if the only non-zero rates are those that appear 
in (3.28) and F^^, (a,/3 G {0,1}). In order to find the degeneracy of the 
ground state, notice that if we set 

1 -poi r02 ri2 

r^f = i> ^ = ^ = ^ = .^ (3.29) 

1 Q,^ 01 ^ 02 12 

we obtain the (2, 1) quotient of the Hecke algebra |1C], and by Schur-Weyl 
duality, the Hamiltonian is Uq{su{2\l)) symmetric. The ground state wave- 
functions are proportional to the monomials symmetrized using the relations 
(3.28,29) and are therefore given by the representation of the superalgebra 



which corresponds to a Young diagram with L boxes in one row |34]. 



A physically important example is the case in which we consider hopping 

11 

Q/3 



and interchange of particles; we set the rates to zero in the previous 



example to get the three relations 

^10 p21 p02 

I^o^i = ^^1^0, ^1^2 = T^D2Di, D2D0 = :^DoD2. (3.30) 
-'- 10 ■'-21 -'-02 
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We have checked whether the Hamiltonian thus defined satisfies the neces- 
sary condition as prescribled by Reshetikhin (eq. (3.9)) for being integrable. 
We found that in order to be so, the rates must satisfy the relations (3.29) 
and we get the (3, 0) quotient of the Hecke algebra and the Hamiltonian 
is Uq{su{3)) symmetric. This result is important because it shows that al- 
though the ground state of the Hamiltonian can be obtained by algebraic 
means, the spectrum of the Hamiltonian cannot, unless supplementary con- 
ditions are satisfied. We have also checked other cases with the same result. 
If the relations (3.29) are satisfied, the ground states again give the repre- 
sentation of Uq{su{3)) which correspond to the Young diagram with L boxes 
in one row [^ ]. 

We can also consider the case where the only rates that are non-zero are 

10 5 01' '-20 ^'l'^ 02- 



r^o, Fq]*, and rg2. The only two relations that survive are 



T%DoDi = T'olDiDo and r^^Z?o^2 = r^2^2l?o, (3.31) 

and the degeneracies are huge since any sequence of the state indices 1 and 
2 is invariant under the dynamics and can be computed using the results of 



ref.m 



Finally, one can have only one relation 

DiDj = 0, (3.32) 

if the only non-zero rates are T^-'^ for any a, [5. 

Since we have only looked at particular cases, it is most likely that we 
have not exhausted all possible quadratic algebras arising from 9x9 intensity 
matrices. We have not been able to find other examples which give non- 
trivial (non-scalar) representations for the relations emerging from a choice 
of intensity matrix. The rules of this game are however very simple and 
can be applied by the reader to find other examples. In the 9x9 intensity 
matrix, one takes several diagonal blocks of smaller intensity matrices with 
non-zero rates and columns of possibly non-zero rates of the form depicted 
below. 
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Notice that there are no rectangular boxes which span entire rows of the 
matrix T above. 

This intensity matrix T gives quadratic relations which then have to be 
checked for consistency for the 3-site case. We can do the same with another 
intensity matrix V := P^yVP^y, where P^^ is a permutation matrix, 

where tt^j.,!^^. are defined before eq.(1.15). The same procedure applies for 
any number of species N . 

To sum up, we have shown that there are a number of quadratic algebras 
which allow the computation of the ground state wave-function. What is 
missing is a complete classification of the algebras. What is also missing is an 
understanding of the origin of the degeneracies. Are the ground state wave- 
functions some representations of some algebra? In the special case where 
the Hamiltonian densities are generators of a Hecke algebra, the answer is 
known, but not in the general case. 

We would also like to emphasize that the ground state wave-functions 
thus computed are non-trivial in the sense that connected correlations func- 
tions in these states are non-zero. One simple example is already known. 
This is given in eq. (3.18) [^^. The method used in [36| to compute the 
correlation functions was different. Other cases can be considered using the 
results of this section. 
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4 Linear chain with left end open and right end 



closed 

As discussed in section 1, the ground state wave-function is given by the 
expression 

L / N~l \ 

\ps)) = (oi n E D^.^^. ' (4-1) 

k=l \Mfc=0 / 

where the D^'s satisfy one of the polynomial algebras described in the last 
section and a new condition (see eq. (1.22)) 

N-l 

Y,Li{Q\Dp = G, (4.2) 

/3=0 

where is an intensity matrix with cofactors C13. Wc consider the cases 
N = 2 and N = 3 only. The following conditions on the D matrices are 
obtained depending on the structure of the boundary intensity matrix. In 
each case, only the non-zero rates are written down. 
N = 2 



1. 

{0\D^ = ^{0\Do (4,L?^0) (4.3) 

2. 

{0\D, = (Lj / 0) (4.4) 

3. 

(0|i^o = {Ll^ 0) (4.5) 

N = 3 

For the case where at least one principal cofactor, e. g. jCq, is non-zero 

(0|i^„ = ^(0|i^o (a = 1,2) (4.6) 

and one obtains two constraints. If all the principal cofactors are 0, only 
two rates are non-zero. We can distinguish the following cases: 
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1. 



2. 



3. 



4. 



5. 



6. 



{0\D, = ^{0\Do 4,L$^0; (4.7) 



(0|L>2 = ^(Oli^o LlL^^^O; (4.8) 
{0\D2 = ^{0\Di LiLj^O; (4.9) 

{0\Di =0 (Lj + L^) 7^ 0; (4.10) 

{0\Do = (L? + L") 7^ 0; (4.11) 

(0|Z^2 = (lI + lD^O. (4.12) 



Notice that for the last six cases, one has only one constraint. The con- 
ditions (4.3) to (4.12) have to be compatible with the polynomial algebras. 
In this way one can, in some cases, find relations between the bulk rates 
only. In most cases the degeneracy of the ground state is lifted. 

The marriage of the bulk algebra with the boundary conditions (4.3- 
12) is a straightforward mathematical exercise. Therefore, we shall confine 
ourselves to models of asymmetric diffusion where the only non-zero rates are 
and give examples of how the concentration profiles can be computed. 
These arc the cases for which the relevant algebras are given in eq. (3.18) 
for iV = 2 and eq. (3.30) for = 3. 

We shall start with the N = 2 case. We have the algebra 

DqDi = q^DiDo (4.13) 

and one of the conditions (4.3—5). The conditions eq. (4.4) and (4.5) give 
a trivial ground state: (O|-Do i^'i- ^■^) i^l^i i^l- 4-5). We shall rewrite 
eq. (4.3) as follows: 

(0|L>i = /x(0|Do. (4.14) 

There are no constraints connecting and Assume that Di indicates 
the presence of a particle on any particular site and Dq the lack of one, i.e. 
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holes. The concentration of particles at a distance k from the left boundary 
for a chain of length L is Q 

c{k) = ^(0|C7'=-^L'lC7^-^ (4.15) 

where 

Zl = (0|C7^, C = L»o + ^i. (4.16) 
A straightforward calculation gives 

= 1 + 42(1-.) ■ (4-17) 



c(A:, /) = ^(0| (&-^D^&-^-^D^C^-^\ - c{k)c{l) (4.18) 



The connected two-point correlation function 
1 

vanishes. 

An interesting case occurs for the algebra 

DiDq = (4.19) 

which corresponds to T^q = and consider the boundary condition = 
(eq. 4.4). Upon inserting these conditions into eq. (4.15) we obtain 

(0| Ej-.i Oijof-' 

and we cannot go any further because our algebraic rules do not allow us to 
evaluate this expression. We leave it as an exercise to the reader to find the 
physical reasons why the method appears not to work in this case, and to 
do the calculation properly. 

We now consider the = 3 states problem and illustrate the method in 
the case of eq. (3.30), which we rewrite for convenience: 

qoDiD2 = D2Di. qiD2Do = DoD2, ga^oA = I?iA, (4.21) 

where 

pl2 7^20 rOl 

10 - 7^21' 91 - ^^02' 92 - T^- (4.22j 
■■-12 ■'-20 J- 01 
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We consider the boundary conditions (4.6) and assume all the rates in (4.22) 
to be non-zero. Setting := Ca/ Cq, we have 

(0|Di =/ii(0|i?o and (0|Z)2 = M2(0|I)o. (4.23) 

As we shall now show, consistency conditions on the rates show up. Let 
us first compute 

go(0|DiL'2 = gom(0|-DoD2 = go9im(0|-D2-Do = gogi/^iA*2(0|£'^; 
{Q\D2Di = fi2{0\DoDi = {fi2/q2){0\DiDo = (Mi/i2/'72)(0|L'g. 

(4.24) 

From the first equation in (4.21) and equations (4.24) we get 

qoqiq2 = 1- (4.25) 

Note that this condition constrains only the bulk rates. The concentra- 
tion profiles for particles "1" and "2" ("0" denotes vacancies) can easily be 
computed ||5| as 

ci{k) = ^(0|C'=~iZ?iC7^-^ C2{k) = ^(0|C'=-il)2C^-' 

Zl = {0\C^; C = Do + Di + D2. (4.26) 
Using equations (4.21), (4.23) and (4.25) we obtain 

= ^ = -«r'- ("-27) 

1 + /^1^2 +tJ'2qf C2{k) /X2 

We are not aware of the existence of another method which would give us 
the expressions (4.27). The connected two-point correlation function is zero. 



5 The open chain three-state exclusion models 

This is a much more difficult problem. We limit our investigation to exclu- 
sion models in which the only non-vanishing rates in the bulk are 

9»i3 ■■= r^f (5.1) 

and make the assumption that the matrices X^'s in equation (1.12) are c- 
numbers, which we shall call x^. As we shall see below, the choice = 



24 



is possible for pathological choices of boundary conditions only. The choice 
(5.1) is not only justified by the interest in these physical processes but as 
shown later, (see the discussion after eq. (5.24)) in this way one obtains 
quadratic algebras which are understood |l^ . 

We are looking therefore at Fock representations of the algebra 

501^0^1 - gioDiDo = xqDi - xiDq 

912D1D2 - g2iD2Di = X1D2-X2D1 (5.2) 

920D2D0 - = X2D0 - X0D2 

with boundary conditions on the states (0| and |0) 

{0\{x, - Li^D^,) = and (x, + i?(;Z)^)|0) = 0. (5.3) 

As shown in ref. Q, for N = 2, one has Fock representations of the algebra 
for any bulk and boundary rates. The situation is going to be very different 
for = 3. The reason is that the number of boundary conditions deter- 
mined by the matrices and is in general too large. For each of L and 
R one has to consider all the cases enumerated in the previous section (see 
eqs. (4.6—12)). We now make some transformations which exploit the fact 
that the x^'s are still free parameters. 
Let us first consider words of length 1: 

{0\D,\0) = 5^. (5.4) 

From equation (5.3), it follows that 

{L^ + R^)6^ = (5.5) 

and 



1 

2 

We have to keep in mind that 



X, = -{Li^ - R^)5^. (5.6) 



s;: = l;; + r; (5.7) 

is an intensity matrix with cofactors This implies (using eqs. (2.3-6)) 
that 

= (5.8) 

and ^ 

= '2^^'^ ~ '^1')^^' (5.9) 
= aiii'R^-R^.C^). 
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Here ^ is an arbitrary parameter and 1Z^{C^) are the cofactors of the in- 
tensity matrix R^{L'^). The expressions for the various cofactors can be 
obtained using eq. (2.7). Since and are intensity matrices, it follows 
from eq.(5.3) that 

2 

^x^ = 0. (5.10) 

It is interesting to note that if we have a Fock representation, the ratio 
of the currents of particles 1 and 2 is already known. It is easy to show that 



where 



This implies 



C=Y.D^.. (5.12) 

At=0 



J2 X2 



where xi/x2 is given by the boundary conditions (see eq. (5.9)). 

In deriving (5.11) we have used the following expression for the current 
density §, | 



ji = qioDiDq + gi2DiD2 - goiD^Di - g2iD2Di 
= -{xo + X2)Di + xi{Do + D2) = xiC, 



(5.14) 



and a similar equation for j2. We also used the definition 

J^ = in). (5.15) 

It is easy to find the one-dimensional representations of the algebra (5.2). 
Using eq. (5.4), we get the following conditions on the rates: 

501 - 510 = yo - 2/1 

512-521 = ^1-52 (5.16) 
520 - 502 = 2/2 - 50 



where 



5m = / (5.17) 
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are given by the boundary conditions (see eqs. (5.8) and (5.9)). Notice that 
eqs. (5.16) give one relation between the bulk rates: 



501 + 512 + 520 = 510 + 521 + 502 (5.18) 

and two relations between boundary and bulk rates. 

Before we start looking for the cases in which the Fock representations 
of the algebra (5.2) exist, it is useful to bring it into a different form. We 
denote 

2/1-52,-1 52-5o,T 50-51,, 

Uo = hl-Qo, ui = \-l-qi, U2 = hl-g2 (5.19) 

521 502 510 

and also 

vio = l-q2 + (50/510) VQi = l-q2 - (51/510) 

U21 = 1 - go + (51/521) t^i2 = 1 - go - (52/521) (5.20) 

^^02 = 1 - gi + (52/502) V2Q = l- qi- (50/502)- 

The g^'s and y^s in the above are defined in equations (4.22) and (5.17) 
respectively. We also define new generators //^ 

D^ = 5^{l + H^). (5.21) 

With these new definitions, the algebra relations and the action of the 
generators on the vacuum are 

q2HoHi - HiHq = U2 + vioHi + -Uoi-f^o 

qoHiH2-H2Hi = uo + V21H2 + V12H1 (5.22) 
qiH2Ho - H0H2 = ui+ V02H0 + V20H2 

with 

(0|Fo|0) = (0|Fi|0) = {0\H2\0) = 0, (5.23) 



iCo/5omHi = (£i/5i)(0|/7o; iCo/5o){0\H2 = (€2/52)^ 
(7^o/<5o)Fl|0) = (7^l/^l)Fo|0); {no/So)H2\0) = in2/52)Ho\0) 



(5.24) 

The algebra (5.22) is of a special case of the ones defined in eq. (1.24) 



[12, 13 1 . For the latter one can show that For the latter one can show that if 
the coefficients are generic, then for > 3 the algebra is finite dimensional. 
For example there are 28 independent monomials for = 3. Therefore one 
can already understand that the = 3 problem is different from the one for 
N = 2. The next point to note is the role of boundary rates for = 3. They 
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appear explicitly in eq. (5.24) but also through the y^'s in the coefficients 
of the algebra (5.22) (see eq. (5.20)). 

In the next section we ask the question if the algebra (5.22) has represen- 
tations in the generic case (when the principal cofactors of the boundary ma- 
trices are non-zero) aside from the one-dimensional one, which corresponds 
to -u^ = and = 0. We shall show that one can have representations 
of dimension at most 2. In general, it is not possible to have higher dimen- 
sional representations because the number of constraints coming from the 
algebra relations is larger that the number of boundary and bulk rates. In 
the subsequent sections, we will make a systematic investigation of the cases 
where some or all of the principal minors vanish. 



6 Generic 3-state exclusion model: boundary in- 
tensity matrices have non-zero minors 

In this section we consider the question when the algebra (5.22) has Fock 
representations defined by eqs (5.23—4) under the assumption that the Tl^i 
and are non-zero. This is a very technical section and the reader not 
interested in the method of answering the question can proceed directly to 
the end of the section where the result is given. For our specific purpose 
it is convenient to rewrite the algebra in a slightly different way using new 
notation. 
Let 

\ ^0 T^i So 2 ( n 1 o ■ 1 o^ 

qv = Qj,; \ = -r^^ m» = ^:r' «i = — > (^^ = o, i,2; ? = 1,2) 

^0 Oi /■Co Oi Hi 

(6.1) 

and instead of the generators Hv use defined as follows: 

Do = 5o{l + Lo), Di = 5^{l + X^L^), = 52{l + H2L2) ■ (6.2) 

The algebra relations are 

QILqLi - LiLq = 1712 + nioLi + noiLo 

QIL1L2 - L2L1 = mo + n2iL2 + ni2Li (6.3) 

Q1L2L0 - L0L2 = mi+ 7102^0 + 7120^2 
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and the conditions for a Fock representation are 

(0|Li = (0|Lo, (0|L2 = af(0|Lo, 
L2IO) = Lo|0), Li|0) = aiLolO), (6.4) 
(0|L,|0) = 0, (z^ = 0,l,2) 

This algebra depends on 14 parameters, and after setting the time scale, we 
are effectively left with 13. For the physical problem, the m,/'s and n^i^'s 
are related to the rates: 



mo = 










mi = 


/"2 V 


m2 = 





V 512 J 



^^ + 1-Qn> (6-5) 

502 

gio 



and 




nio = l-Q2-\ > '^oi = 

1 / \ 
no2 = — (1-Q? + ^), n2o = (1-Qi-^], (6.6) 

M2 V 502/ 

1 / 1 r,2 . yi\ 

n2i = — 1 - Qo -I ' ^^12 = 

Ai V 521/ M2 

At this point two strategics arc possible. One can look for vacuum ex- 
pectation values of words of different lengths and find consistency conditions 
for the 14 paramenters; or, one can look for matrix representations of the 
algebra and in this way obtain consistency conditions among the parame- 
ters. We shall do both - we shall take matrix representations and indicate 
which of the consistency conditions on the matrix elements also come from 
specific vacuum expectation values of words of given length. In order to do 
so, it is convenient to work with paramaters Pi^^, and introduced below 
instead of the m^ and n^i^ introduced earlier. 

mo = <5oaia2-^((5o/(ai«2))ro, 

mi = QiaiN{Qiai)ri, 

m2 = Q2C(2NiQ20(2)r2, 

noi = Q2a2N(Qja2)poi, 

mo = Q2N{Qla2)pio, (6.7) 

n2o = QiN{Qlai)p2o, 

no2 = QiaiN{Qlai)po2, 

n2i = Qooi2N{Ql/{aia2))p2i, 

ni2 = QooiiN{Ql/{aia2))pi2- 
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The constraints on the parameters we shall obtain in what follows look much 
simpler in terms of the new parameters. In eqs. (6.7) and in the what follows 
we use the notation 



1 



M{x) = x + - and N{x) 



1 

X 



{61 



We consider a three-dimensional representation of the algebra eq. (6.3). We 
make this choice since we are going to prove that in the the generic case, 
the two-dimensional representation is the largest we can have. Using the 
boundary conditions (6.4) we can bring the matrix representation of the 
generators to a tridiagonal form by means of a similarity transformation 
by an orthogonal matrix, just as it was done for the N = 2 case in One 
can show that with 

(0| = (1,0,0) and |0) = (1,0,0)^ 
(where T denotes transposition) the generators then take on the expressions 





/ 










Lo = 












I 





vr2 


^2 



/ 



Li 



a 



Li 



V 



2 rr 






/Ti 





(1) 

a2c\ 











1) 



(6.9) 



/2 



(2) 



(2) 

aiC2 



J 



Notice that a one-dimensional representation is obtained if /i = 0, a 

two-dimensional one if /2 = 0, /i 7^ and a three-dimensional one if /i 7^ 

(7) 

0, /2 7^ 0. fi and c- are still to be determined from the algebra relations, 
the Fock conditions (6.4) having already been taken into account. 

We now insert the 3x3 matrices (6.9) in the algebraic relations (6.3) 
and from the equality of each matrix element (i,i) on the left- and right- 
hand sides, get the following equations. From the (1, 1) matrix element, we 
obtain 

ro = ri = r2 = /i. (6.10) 
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From the matrix elements (1,2), (2, 1) and (2, 2) we get 

cf^ = M{a2Q2)pio + M{Q2)poi = M{Qi)po2 + M{aiQi)p2o 

= M{Q2)pio + M{a2Q2)poi = M(Qo/(aia2))P2i + M(Qo)Pi2 

cf^ = M{aiQi)po2 + M{Qi)p2o = M(Qo/(aia2))pi2 + M(Qo)P2i 

(6.11) 

NtfuL ^ ; - /.+.?o+.^. + M(gi.2).o...o 
Nl!^^li^) f^ = A+^?2+pii + M((a,a2)/Qg)p2m2 (6-12) 
J0^h = /r+pio+P^2 + M(Qf«,).02.20. 

Let us pause here for a moment. Equations (6.10) and (6.11) can also 
be obtained from the vacuum expectation values of words of length two 
and three. The 1-dimensional representation is obtained from the condition 
/i = and the 2-dimensional representation is obtained by setting /2 = 
in (6.12). Prom eqs. (6.11) and (6.12) the p^u are completely determined 
as are the (eq. (6.10)). This means that in the algebra (6.3), the free 
parameters are ai,a2 and (3,^,// = 0, 1,2. The next logical step would be 
to solve the equations and look for solutions with positive rates. This is a 
difficult exercise. For a given physical problem, where some conditions on 
the rates are given, the problem is simpler since there are fewer parameters. 

We now go on with the rest of the matrix elements (1,3), (2,3), (3, 1), 
(3,2) and (3,3). We get: 

QlQiQl = 1, (6.13) 

= i^i^Q4) ((^(«iQi) - NiQl))p2o + M(Qi)iV(Qiai)po2) 

= -^^^ (malQi) - NiQl))p,o + M{Q2)NiQ2a2)poi) 
(1) ^ MiQ2^ 
N{a2Qi 



(Af(a^Q^) - N{Qi))poi + M{Q2)N{Q2a2)pi^ 



- Kr. \\ {^^^^^ - ^(^0))P21 + M{Qo)N{^)p,2) 

c?^ = ((^(«?Qf) - Ar(g?))po2 + M(Qi)iV(Q2a2)p2o) 

= ((iV(^) - ml))P^^ + MiQo)Ni^)p2,) 

jy-/_vo_\ V aia2 aia2 ^ 

^ aia2 ' 

(6.14) 
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and 



= -fl + 7V(^^(^(^l)4'^4'^ - ^(«lQ?)feocf +P024'^)). 

(6.15) 

Equations (6.11) and (6.14) give a system of six homogeneous equations for 
the p^i, with determinant genericahy different from zero. In order to have 
a non-trivial solution for the Pf^u^s, we must set this determinant to zero, 
thus introducing an extra condition on the rates. This leaves one of them 
(say pio) free. The same equations determine c[^^ (i = 1,2, = 0,1,2). 
We are thus left with eight equations: three from (6.12), one from (6.13), 
three from (6.15), and one determinantal condition for the seven unknowns 
f2,Pio^ Qii and Oj. (We can take /i = 1 in eq. (6.10) to set the time scale.) 
This implies that we cannot, generically have a 3-dimensional representation. 
This also implies that the algebra (6.3) with the Fock conditions (6.4) can 
have representations of dimensions 1 and 2 only. The conditions on the 
parameters are not neat but this was so even for the case = 2 [^. 

We now come back to the physical problem. The algebra (5.22) contains 
9 parameters - qo,qi,q2 and the ratios (yoMo), {vo/gm), (yi/gio), (yi/521), 
(^2/502) and (2/2/521); the boundary conditions (5.24) depend on 4 others, 
the Aj's and /ij's for i = 1, 2. This makes a total of 13 parameters. Therefore 
for the physical problem, the case of 2-dimensional representations stays 
interesting since the algebraic varieties on which the solution set exists is 
not trivial. 



7 The two non-zero cofactors case 

In the last section we considered the generic case where all the cofactors 
and TZ^ are non-zero and we found that the algebra (5.22) had Fock 
representations defined by eqs. (5.23) and (5.24) of dimension 1 and 2 only. 
We should have continued our study and looked at all cases when some or 
all of the cofactors vanish. In this section we shall confine ourselves only to 
the case 

£0 = ^:2 = ^0 = 711 = 0, £l/0,7^2/0. (7.1) 
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We had in mind the physical apphcation of ref. The interested reader 
can easily duplicate the calculations presented in this section for the other 
cases. With (7.1), the Fock representations have to satisfy the following 
simple simple conditions 

(0|Fo = (01^2 = 0; Fo|0) = Fi|0) =0, (0|F,JO)=0. (7.2) 

In order to find out in which cases the algebra (5.22) has Fock representa- 
tions defined by eq. (7.2), we take vacuum expectation values for monomials 
of different degree. 

Monomials of degree 2 ( consider the algebra relations between (0| and |0)) 
give: 

ui = U2 = 0. (7.3) 

Monomials of degree 3 (consider {0\HiHqH2\0) and commute Hq through 
to the left and the right) give 

vio = V20 = V. (7.4) 

Monomials of degree 4, {0\HiHiHoH2\0) and {0\HiH2HoH2\0) give 

qovoi = q2V2i and qoVo2 = QiVu (7.5) 

respectively. We have to distinguish several cases 
• a) vi2 and V21 are non-zero. 

Monomials of degree 5 (consider {0\HiHqHiH2H2\0)) give 

91 = ^2 = q- (7.6) 
Taking into account relations (7.3—6) the algebra (5.22) becomes 

qHoHi-HiHo = vHi + ^V2iHo 
a ^° 

qH2Ho-HoH2 = —V12H0 + VH2 (7-7) 

qo 

qoHiH2 - H2H1 = uq + V12H1 + V21H2 
It is useful to change variables 

-Ho = -vN , Hi = — A, H2 = — o, Z2 = (7.8) 

Vl2 V21 1'121'21 
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in terms of which, the algebra becomes 



a) 22 (aB - —BA 

\ qo 

h) AN - qNA 



l+A + B 



A-^M 
B-^M. 

Z2 



(7.9) 



c) MB - qBM 

This algebra has the representations 

(0|AA = = AA|0), (0|e = = ^|0). 



(7.10) 



We can now use the results of ref. ^ to find the following representation 
of (7.9a): 



A 



where 



/ ai /i 0... \ 

aa /2 0... 

ag /s... 

04... 



B = A^ 



(ax 0... \ 

/i aa 0... 

h a-s 0... 

/3 04... 



a„ = Z2 ^{n - 1}a, 
1 — A" 

= " ^ 1' iO}A = 1 



f^ = z^^{n}x{z2 + {n-l}x) 
A = 



(7.11) 



(7.12) 



We have introduced the symbol Z2 in keeping with the notation of ref. 
[^. A'^ is the transpose of matrix A. There are two cases when we can have 
representations for J\f in eqs. (7.9): 



a.l) 



when 



qo, Vol = V21, V02 = V12 



f hi . . . \ 
/l2 ... 

/13 ... 



with 



hn = {n- 1}^ 



(7.13) 



(7.14) 



(7.15) 
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and 
with 



a.2) q = 0, go / 0, voi = V02 = 



hn = {n- l}q^O = 1 - <5nl 

• b) We consider the case (see eq. (7.5)) 

Vol = V02 = Vl2 = V2I = 0, Vio = V20 = V, Ul = U2 = 0. 

Monomials of degree 5 give (see eq. (7.6)) 

q\ = q2 = q- 

With (7.18) and (7.19) the algebra (5.22) becomes 



qHoHi - 


HiHq 


= vHi 


qH2Ho - 


H0H2 


= VH2 


qoHiH2 - 


- H2HI 


= Uq. 



We make the transformation 
Ho = -vAf, Hi 
and get the algebra relations 



'^A, H2 = J^B, X = - 



qo 



qo 



AB - XBA 
AN-qMA 



1 

A 
B 



and 



NB - qBM 

A\0) = M\0) = = {0\B = {0\M. 



(7.16) 
(7.17) 

(7.18) 
(7.19) 

(7.20) 



(7.21) 



(7.22) 



(7.23) 



The Fock representation of this algebra (matrices satisfying (7.22) and (7.23)) 
is: 



A = B'^ 



/ 51 . . . \ 
52 




V ■■ 



(hi ... \ 

/l2 
/t3 



■• / 



V ■■ 



where 



9I = hn = {n- 1} 



(7.24) 



(7.25) 
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• c) We consider the case (see (eq.(7.5)) 

V21 = vqi = 0. 
The algebra (5.22) again with qi = q2 = q becomes 



qHoHi ~ HiHq 



vHi 

q 



qH2Ho — H0H2 = —V12H0 + VH2 

<1Q 

qoHiH2 - H2H1 = ■uo + V12H1 
We make the substitutions 

Ho = -vAf, Hi = ^A, H2 = —B, A 
V12 qo 

to get the algebra 

AB - XBA = 1 + A 

AM~qMA = A 
MB-qBM = B-qM. 

This algebra has the representations 

/0/iO...\ /61OO 
/2 




A 



\ 



B 



/i 62 

/2 h 



\ ■ 



where 
and if 
then 



fn = Wa, bn = {n- 1}; 



c.l) q = qQ = X \ 



where 
If 



/ hi 





. 


. \ 





h2 













hs 




V '■ 






• / 


hk = 


{k- 


1}a. 




C.2) 


Q 


= 0, 
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which imphes V02 = one has to take 

hk = {k - Ijx^o = i - Ski. (7.33) 

In order to help the reader, we shall summarize the results. The algebra 
(5.22) has Fock representations defined by 

(0|i/2 = (0|i/o = = i?o|0) =//i|0) (7.34) 

if {necessary condition) 

ui = U2 = 0, vio = V20 = v, qi = q2 = q, (7.35) 

in which case the algebra relations become 

qHoHi - HiHq = vHi + vqiHq 

qH2Ho-HoH2 = V02H0 + VH2 (7.36) 
qoHiH2 — H2H1 = uo + V12H1 + V21H2 

This algebra has Fock representations in five cases 
1. 

vqi = V21, vq2 = V12, qo = q- (7.37) 

With the notations 

Hq = -vN, Hi = — A, H2 = — o, Z2 



V12 V21 ■yi2'y2i 

A and B are given by eq. (7.11-12) with A = q^^ and M by eqs. (7.14) 
with (7.15). 

2. 

vm = vq2 = q = 0. (7.38) 
Same as case 1) except that M is given by eq. (7.14) with (7.17). 

3. 

Vol = V02 = V12 = V21 = 0. (7.39) 

Ho = -vAf, Hi = .f^A, H2 = .f^B, A = -. (7.40) 

V 90 y qo qo 

A, B and J\f are given by eq. (7.24-5). 
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4. ^ 

Vol = ■V21 = 0, q = qo = -. (7.41) 

Hq = -vM, Hi = ^A, H2 = —B, A = - 
^, B and A/" are given by eq. (7.29) and M by eq. (7.31) with (7.32). 

5. 

vo2 = vqi =V2i = 0, q = 0. 
Everything as in case 4) except Af is given by eq. (7.31) with eq.(7.33). 

This closes the problem of Fock representations of the algebra (7.7). We 
now turn to the physical problem. This implies using the definitions (5.19) 
and (5.20) of the parameters and v^i,. We first have to give the boundary 
rates which give eq. (7.1). Using (2.7) we have the solution 

L2 = Lq = i?^ = Rq = 

£1 = L\Ll + LlL\ + LlL% (7.42) 

We now compute and using eqs. (5.8) and (5.9): 
<^o = ^{RqL'i + L'^Rq + -^0-^2) 

(5i = iiR'iLl + R^Ll + R'iLl + Ci) (7.43) 
52 = ^{RlLl + RlLl + R\Ll + n2) 

and 

xo = i{LllZ2-RlCi) 

XI = i{L\n2 + {Rl + R\)Ci) (7.44) 
X2 = -^{{Ll + Ll)n2 + R\Ci). 

We now implement the conditions (7.35) to get 

yo = 0, 2/1=510(1-9), y2 = -502(1 - 9)- (7.45) 

We have to consider the 5 cases separately. 
1. Prom (7.37) and (5.20) we get 

yi = = y2. (7.46) 

Taking into account the positivity of the rates, the definition (5.17) 
and the relation (5.6) this can be excluded. 
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2. From eq. (7.38) and (5.19) we have 

520 = = 501, 2/0 = 0, yi = giQ, y2 = -go2- (7.47) 
Using the definitions (5.19) and (5.20) we get 

^^10 = V20 = 1 

and 

, ,502 ^ , 5io -, , 5io + 502 

vi2 = l- qo-\ , V2i = l-qo-\ , uq = 1 - qo -\ , 

521 521 521 

and the following constraints on the rates coming from the definition 
(5.17) of and from from eq. (7.47): 

^0^2 = Rq^i 



510 



RlLl + RO.Ll + R^^q+Ci (7.48) 
{Ll + Lf )7^2 + RlCi 

3. As a consequence of eqs.(7.39) and (7.45) we find: 

510 = 502, 501 = 520, 512 - 521 = 5io - 5oi, (7.49) 

yo = 0, yi = -y2 = 510 - 501 (7.50) 

and 

'yi2 = = t;2i, UQ = qQ-l, viQ = V2Q = 1 - q. (7.51) 
Prom eq. (7.50) we find the following constraints on the rates: 

LqT12 = RqJ^i 

LIIZ2 + {Rl + R\)Ci {Ll + Ll)lZ2 + R\Ci 



R^Ll + R°Ll + R^Ll + £1 RlL^ + RlL% + RlL\ + 7^2 

= 510-501- 

(7.52) 

4. No solutions if we want to maintain positivity of the rates. 
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5. We find 

521 + 510 = 512 520 = = 501, (7.53) 
50 = 0, yi = 510, 52 = -502 (7.54) 

and 

, 502-510 502 1^ v^r^\ 

VlQ = V2Q = 1, V12 = , Uq = . (7.55) 

521 521 

The conditions on the rates can be derived from eqs. (7.54) and (5.17). 
Notice that this case is a speciahzation of case 3. The representation 
is different because Z2 defined in eq. (7.8) and used for case 3 diverges 

{V2l = 0). 

In Appendix A we consider the case of a CP-invariant stochastic process 
[^. In this problem, the index "0" denotes vacancies, "1" is the index for 
particles and "2" for anti-particles. This implies supplementary conditions 
on the rates. If we are interested in the case where the CP-symmetry is 



explicitly broken 38 1 one can use the solutions corresponding to the 
cases 2), 3) and 5) described above. In all three cases the D^s have the 
form 

Da, = XaEa {Xa / 0), Da = Ea {Xa = 0) (7.56) 

where E^, which can be extracted from the formulae given above and are 
given explicitly in ref. |2^, depend only on the bulk rates. The boundary 
rates enter in the expressions of Xa and through the fact that they have to 
"match" the bulk rates (see eqs. (7.48) and (7.52)). This observation allows 
us to see the important parameters of the problem. 

Before closing this section, let us note that the case where only the 
cofactors Cq and TZq vanish can be treated in a similar way. In this case, the 
algebras satisfied by the A and B matrices (see equations (7.9a) and (7.28)) 
contain A"^ and terms as well. Fock representations for this case are also 
known B. 



8 Lower rank boundary intensity matrices 

In sections 6 and 7 the Fock representations were defined by four conditions 
(see eqs. (5.24), (6.4) or (7.2)). The number of conditions depends on the 
rank of the boundary intensity matrices. In principle, one should examine 
all the possibilities enumerated in eqs. (4.6—12) for both left and right 
intensity matrices. This is a long exercise. The main point is that once 
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either LJ^ or or both have all the principal minors to be zero, the number 
of conditions is smaller and one has more freedom for the bulk rates. We 
shall give only four examples, 
(a). L\ 7^ 0, all other cofactors zero. 

We take the following to be the boundary intensity matrices: 

(-L\-L\ Ll \ l-R\ Rl 

L = Li Lj \ and R=\ Rl -R^ 

\, -Lj-Ll I \, 

We set ^ = 1 in eqs. (5.8) and (5.9) to get 

xq = —RqCi, 
x\ = RqCi, 



(8.1) 



X2 = 

6o 



rI{lI + lI), ^ 

Si = R^iLl + LD+Ci, 
62 = R^Ll 



We get several constraints on the various parameters vij and Ufc that 
appear in the algebra (5.22) and are given in terms of the bulk and boundary 
rates in (5.19-20). From the condition X2 = 0, and by sandwiching (5.22) 
between (0| and |0), we obtain 

V21 = no; V2Q = ui = 0. 

Prom {0\HiH2HoH2\0) we obtain 

■vio = 

from which, on calculating {0\HiHoH2\0) , we get qi = go = q- The correlator 
{0\HiHoHiH2\0) gives 

92 = 1 + — (9- 1), 
uo 

so we end up with 

Xq 

T- = (501 - 510) = (502 - 520) = (521 - 512) ,^ 

c»o (8.3) 
V12 = V02 = (1 - 9) 

The algebra is of a similar form as case c) of section 7 (eq. 7.26), apart 
from a relabelling of the indices and the presence of an additional constant 
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term: 

qHiH2-H2Hi = V21H2 + {1 - q)Hi + V21 

qH2Ho-HoH2 = il-q)Ho (8.4) 

q2HoHi - HiHq = vqiHq + U2 

Introducing 

Hq = —A, Hi = V21B, H2 = {q- l)M, A = and 7 = - — (8.5) 

Vol goi 

we get (compare with eq. (7.28)): 



AB - XBA = 7(1 + A) 
AN-qMA = A 

MB-qBM = B-M+ ^ 



(8.6) 



l-q 

Also, instead of four Fock conditions as in section 7, we have only three: 

(0|^ = {<d\M = 0, B\<d) = a^|0), (8.7) 

with 

We have not looked for explicit represenations of this algebra, 
(b) All the cof actors are zero. 

We consider the case when only the following four boundary rates do 
not vanish: 

L'^2,Ll,Rl and R^. (8.9) 

We find from (5.9) 

x^ = (// = 0,1,2) (8.10) 
and the algebra is a known quadratic algebra (see section 3) 

qoDiD2 = D2Di, qiD2Do = DoD2, q2DoDi = DiDq, (8.11) 

with only two conditions: 

dO rO 

Di\0) = -^DolO) and (0|L»2 = ^{0\Do. (8.12) 

All vacuum expectation values of degree L can be expressed in terms of 
(0|Z)o|0)- The appearance of the polynomial algebras in this new context is 
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very interesting since, as shown in section 3, we know plenty of them, and 
not just for the simple exclusion processes. This opens up the possibility of 
a new class of solutions for more general processes with open ends. 

(c) All the cofactors are zero, but the non-zero boundary rates are: 

L\,L%rI and R\. (8.13) 

In this case the algebra is also given by (8.11) and = 0, (// = 0, 1, 2), but 
the generators satisfy the conditions 

c-o 

(0|L>o = and r>i|0) = ^i^olO). (8.14) 

It is interesting to observe that no boundary condition appears for the op- 
erator D2. Using (8.11) and (8.14) it is simple to see that the ground state 
will have a the single word I?2 , which corresponds to the lattice filled by 
particle "2". This can be expected physically from (8.13) since the species 
"2" is never removed from the chain at either end, unlike the species "0" 
and "1", and is created at the left end with a certain rate L^^Q. Therefore, 
when the time goes to infinity the system will end up with only particles of 
species "2". 

(d) All cofactors are zero, but the non-zero boundary rates are: 

L\,lI,rI and R^. (8.15) 
The algebra is again given by (8.11) but the conditions are 

(0|Z)o = Dq\0) = 0, (8.16) 

and no conditions for Di and D2. In this case we have L degenerate ground 
states expressed in terms of the words DiDJ^^"" (n = 1, 2, . . . , L) . This is ex- 
pected physically from (8.15) since in the large time limit, we should expect 
no "0" particles and the particles "1" and "2" are conserved separately. 



9 Conclusions 

This was a long journey. We will first list the mathematical results which are 
independent of the physical applications. In section 3 we have given a list 
of polynomial algebras defined by a set of homogenous quadratic relations 
(see for example eqs.(3.25-27)). In sections 6 and 7 we have given Fock 
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representations for the algebra defined by eq.(5.22). Those are a special 
class of the inhomogenous quadratic algebras defined by eq.(1.24). The 
representations are one dimensional (see eq.(5.16)), two-dimensional for the 
"generic" case described in section 6, or infinite dimensional (see the five 
special cases of the algebra (7.36)). These algebras can be useful for physical 
applications other than for the case of open chains considered in this paper. 
For example we used them ourselves while studying stochastic processes on 
a ring ||2^. As for the physical applications relevant to this paper, we have 
looked for steady states of stochastic processes defined by bulk and boundary 
intensity matrices (see eq. (1.1)). The probability distribution of the steady 
states is given by the expression (1.13) for chains with open ends. 

In our approach the same expression without |0) is used for a left open 
end and a closed right end. In the case of closed chains we have to drop 
the (0| and |0) symbols. Polynomial algebras are defined by taking the X's 
equal to zero in eq.(l.ll). They exist only for special choices of the bulk 
rates. A given polynomial algebra can be used in three places: for closed 
ends (see section 3), for the left end open and the right end closed when 
supplementary restrictions come from the first equation in (1.12) with the 
X^s are taken to zero (see section 4) and even for the case of a chain with 
open ends (see the example in section 8, eqs.(8.11-12)). As far as we know 
it is the first time that the problem with one or two closed ends is solved 
using algebraic methods. 

Most of this paper is dedicated to the problem of the chain with open 
ends where we have taken c-numbers for the X's in eqs.(l.ll) and (1.12). 
We have confined ourselves to the simpler problem of exclusion processes. 
In solving this problem a crucial point is the rank of the boundary intensity 
matrices since they dictate the number of conditions which define the Fock 
representations. In the 3 state problem, they can be four (see sections 6 and 
7) three or two (see section 8). In the case of four conditions an important 
role is played by the number of cofactors of the boundary intensity matrices 
which vanish. All in all we have given many examples which can be used 
for physical applications. It is obvious how the present method may be 
generalized for four or more states problems. The case of periodic boundary 
conditions is considered in p^ . 

After writing this paper, one can ask where lies the mathematical beauty. 
Probably in the simpler cases where polynomial algebras can be used and 
in the algebras described in section 7. The problem of open boundaries 
with two states depends on 5 parameters; the one with three states on 17 
parameters and there is a price to pay in order to find solutions. In this 
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paper we have never properly exploited the Krebs-Sandow theorem (eqs. 
(1.12-15)) since we have always chosen representations of the X^s to be 
scalars or zero. The reason is that we wanted to have quadratic algebras 
that are understood, like the polynomial algebras or those of type (1.24), 
and we did not find new ones. 

We would like to make one last comment which has to do with the 
connection between finding the ground states by algebraic methods and the 
integrability of the Hamiltonian. As is shown in section 3 several exactly 
integrable Hamiltonians, connected with stochastic dynamics, can have their 
ground state wave function written in terms of words defined by algebraic 
relations. But in general, the existence of an algebraic form of the ground 
state wave function is not a sufficient condition in order to ensure exact 
integrability. 
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A The CP-invariant steady state 



It was recently shown that in steady states CP-invariance can be sponta- 
neously broken [0] and that the phase diagram cannot be obtained by mean 
field methods ||2J]. This is not the right place to describe the physics of the 
problem and the interested reader should have a look at the references. The 
model where this invariance is spontaneously broken is a special case of the 
problem studied in section 7, where the rates we choose are CP invariant. 
Charge conjugation (C) is defined by the operation 1^2,2^1 and — > 0, 
while the parity (P) operator corresponds to interchanging left and right. 
The non- vanishing boundary rates are 

r2_Rl rO_pO rO_pO t2 _ r>l 

^1 — ^2^ — ^2 — ^1^ ^0 — ^0 ( A ^) 

£l=7^2, £o = = 7^o = 7^l = 0. ^ 

The bulk rates are 

510 = 502, 501 = 520, qi = q2 = q, {A.2) 

and gi2 and §21 are arbitrary. 

In ref. a special choice was taken: 

L? = i?^ = Q and lI = rI = /?, {A.3) 

with all other boundary rates set to zero. The non-vanishing bulk rates were 
taken to be 512 = 510 = 502- Here we consider the general case. 

From equations (7.43), (7.44) and (5.17) and choosing = (Lg + 2Lf ) 
we obtain 

5o = Ll 5i = 52 = {L\ + Ll) 

yo = 0, yi = -2/2 = 70X70 + 

+ ^2 

We now use the results of section 7 (equations (7.43)-(7.55)). From eq. 
(7.45) and (A. 4) we have 



tQt2 

L\ + V2 



I r2 



+ ^1 = 510-501- (-4.5) 



We now consider the cases (2), (3) and (5) separately. 
Case (2), 320 = 5oi = 0. 

We use equations (7.47) and (7.48) to get 

7-07-2 

-^1-^0 , t2 



+ Li = 510, (^.6) 
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and 



Vl2 



V21 = 



= 1^20 = 1, 
510 + 521 



512 



521 



Mo 



Vol = V02 = 0, 
_ 2gio + g21 - 512 

521 



(^.7) 



Using eq. (5.21) it is easy to show that 



Introducing 



we have 



[1 + Hi), D2 



[I + H2). 



1 

90 



521 , 510 
and uj = 



512 



512 



{A. 



(^.9) 



V2a; + A - 1 



«1 










U2 


V2 








U3 



/ 1 





V ; ; ; 



{A 



J 

10) 



where 



{^}a 1 + 



{uj + \-lf 



[00 



A 



w 



{2uj - 



1 + A) 
-1}a. 



{n-l]> 



(All) 



1 + A 

Notice that the boundary conditions appear only in (A. 6) and in the 
normahzation factors of Di and D2. They do not change the physics of the 
problem which is governed by lv and A. In other words, taking only the rates 
(A. 3) or the general case (A.l) does not change the physical results. 

In ref. only the case A = was considered. 
Case (3) 510 = 502, 5oi = 520 and gi2 - 521 = 5io - 5oi- 

We make use of eqs. (7.49—52). The bulk rates and the boundary rates 
are related by the constraint: 



510 - 501 



7-07-2 

-^1-^0 



+ Ll 



(A12) 
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We then have 



R 



^0 



1 + R) 



/ ri 

rs 



V : ; : 



Do 



/ si 

S2 

S3 

\ ■-. \ ■-. 



{A.U) 



where 



with 



1 - A*^ and Sk = I - {k - 1}^ 



(^.14) 



(A.15) 



, 520 901 
and g = = . 

512 502 510 

Notice once more that the boundary terms are not essential and in this case 
as well, the conditions (A. 3) capture the whole physics. 
Case (5) (see eqs. (7.53—55)). 

Since 502 = 5i0) we get V12 = (see eq. (7.55)) which together with 
(7.41) brings us back to case (3). 

In conclusion, considering the most general boundary rates (four instead 
of two) does not change the domain of applications of the algebraic approach. 
What we have shown however, is that it can be applied for a larger choice 



of bulk rates ISg] . 
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